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I. INTRODUCTION 



The study of electronic transport in quantum mesoscopic disordered systems has been a topic of interest for a long 
time ^. Due to the random positions of the impurities in such systems, quantum interference effects give rise to 
strong fluctuations in the conductance g from sample to sample. These effects can be observed in a single disordered 
sample by changing, for example, the applied magnetic field |^ since this is similar to a change in the impurity 
configuration of the sample. As a consequence of these fiuctuations a statistical study of the conductance is required. 

Many efforts have been made in order to have a complete statistical description of the conductance in the three 
different regimes of transport ||]: metallic ^ L), where ^ is the localization length and L the typical size of the 
system, insulating ^ L) and crossover L). It is known that in the metallic regime the distribution of the 
conductance P{g) is Gaussian, so the first and second moments, i.e., the average {g) and the variance var((7) are 
enough to describe Pig). It turns out that in the deeply metallic regime vaT{g) is a pure number independent of the 
details of the system known as the universal conductance fiuctuations, depending only on the presence or absence of 
time reversal symmetry and spin-rotational symmetry Q. In the opposite regime of transport (insulating regime), 
the distribution of g is log-normal, which means that In g follows a Gaussian distribution. 

The intermediate regime of transport between the metallic and insulating regimes can be reached, for example, by 
increasing the disorder in a metallic sample. As the disorder is increased the conductance fiuctuations grow in such 
way that va,T(g) becomes of the same order as (g). In this case, the first two moments are no longer enough to describe 
P{g). In fact, the full distribution of g is needed in order to have a good statistical characterization of the electronic 
transport. However, not much is known about the distribution of the conductance in the crossover regime, even in 
the case of a simple geometry like a quasi one-dimensional system or quantum wire (L >> W , where L is the length 
and W is the width of the system), where a smooth transition from the metallic to the insulating regime exists. 

There are a number of numerical simulations in the intermediate regime which show a broad asymmetric distribution 
of g with a flat part for values oi g < 1 and a strong decay for g > 1 |^, || ^, 10| . Also, these numerical results 
have shown an interesting qualitatively similar behavior of P{g) for quasi one, two and three dimensional systems 

With respect to analytical results in the crossover region, the first two moments of the distribution of the conductance 
have been calculated, in fact, for all values of disorder for quasi one dimensional systems, using the supersymmetric 
non-linear sigma model [a model) For the complete distribution P{g) beyond the metal and insulator limits, 

a systematic method has recently been developed by two of us iQ. Using this method, we were able to calculate P{g) 
for quasi one-dimensional systems near the crossover regime, approaching on the insulating side. It was predicted that 
on the insulating side P{g) follows a "one-side" log-normal distribution cut off by a Gaussian for g > 1. Numerical 
calculations have shown a cutoff at 5 = 1 in the crossover and insulating regimes ||, |9[ ^ and a "one-side" log- 
normal distribution for g < 1, in the insulating region [^0[ |ll|. The method reproduces other well known insulating 



'Present address: 23, rue du Loess, F-67037 Strasbourg Cedex (France) 



2 



and metallic limits. For example, using the results of Fig. |l| shows the average and variance of 5 as a function 
of the disorder parameter F = ^/L; the first two moments agree quantitatively well with those from the a model, in 
the metallic and insulating limits. In Fig. ||, P{g) is plotted for a metal and insulator case, using again the results of 
Ip^ , which shows that the two limits are well captured by the method. However, the crossover regime is qualitatively 
correct only on the metallic and insulating sides, suggesting that the approximations made in [ p^ are not as good for 
the complete crossover region. Indeed, numerical simulations have shown that the drop in P{g) at g > 1 appears to 
be exponential ||^, as opposed to the Gaussian cutoff predicted in |Q. It is therefore important to improve the 
approximations made in in order to obtain a better description of the crossover regime. 

In this work we calculate the distribution of the conductance in the crossover regime using the proposed method in 
iQ with some improvements. We show how P{g) changes from a broad highly asymmetric "one-sided" log-normal 
distribution on the insulating side of the crossover regime to a Gaussian-like distribution on the metallic side of 
the crossover regime, as function of the disorder F. The distribution at the crossover regime agrees with numerical 
results, except exactly at g = 1, which may be due to the existence of a singularity ^ which might make our 
approximations less accurate. 



II. SCATTERING APPROACH AND DMPK EQUATION 

In order to study the conductance, we will use the scattering approach to electronic transport, in which the 
(dimensionless) conductance is given in terms of the transmission eigenvalues T„ by 



N 



(1) 



where N is the number of channels or transverse modes of the wire. Due to the random positions of the impurities 
in a disordered sample, the eigenvalues T„ fluctuate from sample to sample and the distribution of g is given in 
general by P{g) — (^S (^g — J2iLi where (...) indicates the average with the joint probability distribution of the 

transmission eigenvalues p({T„}). 

For a mesoscopic quasi one-dimensional quantum wire, the evolution equation for the distribution p(|T„|) as a 
function of the length of the system is given by the Dorokhov-Mello-Pereyra-Kumar (DMPK) equation [Q. For 
systems without time reversal invariance or unitary symmetry, which is the case that we study here, the DMPK 
equation can be written as 



dL 



1 ^ 
N ^ 



n=l 



(2) 



where 



N N 

n n 



(3) 



A„ = (1 — Tn)/Tn and I is the mean free path. The solution of the DMPK equation has been found by Beenakker and 
Rejaei |0 for all degree of disorder and it is given by 



p({a;i}) = CjJ^ ( sinh^ a;^ — sinh^ Xi) JJ^ sinh 2a 



Det 



-k^s/4N 



tanh(7rA:/2)A: 



2m -1 



X P(ife„i)/2(cosh2a;„) 



(4) 



where the Xi's are related with the Ai's by Xi — sinh^ a;^ or with the transmission eigenvalues by Ti = 1/ cosh^Xi; C 
is a normalization factor and P(ik-i)/2 is a Legendre function. 

Using the Fourier representation of the 5-function, the distribution of the conductance is given in general by 



P(9) 



dr 
2^ 



N 



exp 



IT \g 



N 

E 



cosh^ Xi 



(5) 



with p{{xi}) given by Eq. (|^). The calculation of P{g) from Eq. (||) involves a nontrivial N-fold integration. 
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III. FEATURES OF THE CURRENT METHOD 



In the metallic regime, where the transmission eigenvalues T ~ 1 contribute significantly to the conductance, the 
A's {— (1 — T)/T) are very close to each other and a continuum density for A can be assumed. In this approximation 
the universal conductance fluctuations can be derived. On the other hand, in the insulating regime, the A's are much 
larger than the unity (small transmission) and separated exponentially. In this case, the lowest eigenvalue gives the 
most important contribution to the conductance and within this approximation the log-normal distribution for g can 
be derived. But none of these approximations can describe the crossover regime where the eigenvalues are neither 
very close to each other, nor exponentially separated. We will see, however, that it is possible to study P{g) in the 
crossover region by making systematic corrections from the metal and insulating regimes. 

In order to calculate the distribution P{g) in the crossover regime, we propose, following |p^ , to separate out the 
two lowest eigenvalues xi,X2, treating the rest as continuum. In only the lowest eigenvalue xi was separated out. 
We will see that by separating out one additional eigenvalue we can go beyond the region studied in [^. In the 
following we give a brief description of the important features of the current method in order to calculate P{g)- 

The general solution of the DMPK equation (^ can be simplified for the cases I <^ L ^ Nl (metal) and I <^ Nl <^ L 
(insulator). In these cases the integral over k can be calculated analytically and the determinant is given, for example 
in the metallic regime, by a product of the difference of the square eigenvalues The solution for both metallic 
and insulating regimes can be written as: 

p({xj) = |exp[-i/({xj)] (6) 

with 

N N 

i<j i 

The H {{xi}) function can be interpreted as the Hamiltonian of N charges at the positions Xi with a two body 
interaction given by u(xi,Xj) and confinement potential V{xi). Those terms are given by 

u{xi, Xj) = — In I sinh^ xj — sinh^ Xi\ In jx^ — a;^|, (8) 



Irxf — \ ln(a;i sinh 2xi 
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Yxl — i ln(xi sinh 2xi 



■ In: 



(metal) 
(insulator). 



(9) 



The parameter r(= ^/L) is given by F = Nl/L in quasi-one-dimension. F 3> 1 and F ^ 1 correspond to the metallic 
and insulating limits, respectively. F ~ 1 is the crossover regime. We can see fromEq. (§) that the one body potential 
V{xi) in the insulating regime differs by a logarithmic term from the corresponding V{xi) for a metal. In the localized 
regime the transmission is very small, i.e. Xi ^ 1 (since T = 1/cosh^a;) and then the logarithmic term is negligible 
compared to the other terms of V{xi). This means that the solution in the metallic regime might be used for the 
insulating regime as well. Therefore, we will assume that the solution in the metallic regime is also valid for the 
intermediate regime. We remark that it is not always possible to write the general solution p{{xi}), Eq. (^), in the 
form Eqs. (6-7). However, we will check the validity of our approximations by comparing the average and variance of 
g with the respective exact results from the a model for those quantities p^ . 
Using Eqs. (ph and (o), the distribution of the conductance can be written as 



Y roo ^ poo N 
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cosh Xi ] 
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where Z is a normalizing factor. In order to calculate the distribution of g, as mentioned before, we start by separating 
out the two lowest eigenvalues xi and X2 in H(x), Eq. (^: 



3<i<j 



3<i 



with 



Hi.2 = u{xi,X2) + '^u{xi,Xi) + '^u(x2,Xi) + V{xi) + V{x2) 



(11) 
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Now we make the continuum approximation: 



H{xi,X2,X3,a{x)) = V{xi) + V{x2) + u{xi,X2) + / dxa{x)u{xi,x) + / dxa{x)u{x2, x) 



X3 



X3 



2 .1X3 



dx I dx'a{x)u{x,x')(T{x') + / dx(T{x)V{x), 



(13) 



where we have introduced the density of eigenvalues (7{x) which has to be calculated in a self-consistent way and 
subject to the conditions (y{x) > and J^^ a{x)dx = N — 2. 
Defining the "free energy" F as 



F(a;i,a;2,a;3;cr(a;)) = 2H{xi,X2,X3]a{x)) + it 



1 



1 



+ 



a{x) 



cosh^(a;i) cosh^(2;2) 
the distribution P{g) can now be written as a functional integration: 

P(9)^^ TT^'"' dxj dx2 rfa:3 / I?[a(^)]e-^(--^--^^'^(-» 

Joo ^'^ Jo Jxi Jx2 J 



cosh X 



(14) 



(15) 



As in |l^, in order to find the density a{x) self-consistently, we minimize the free energy by taking the functional 
derivative: 5F(xi^X2TX3;a{x))/5(j(x) = 0, which gives the following integral equation for the saddle point density 
asp{x): 



IT 



dx'u{x, x)asp{x, x) — u{x, xi) + u{x, X2) + 2 

cosh X 



V{x). 



(16) 



The above integral equation can be solved when the lower limit goes to zero. In this case, the integral can be 
symmetrically extended to negative values which allows one to invert the kernel and obtain fjgp{x). For the non- 
zero limit, we calculate asp{x) using the following "shift approximation". In terms of the eigenvalues A and A', the 
interaction term in iJ, Eq. (0), can be written as: 



w(A,A') = -i[ui(A,A')+w2(A,A')] 



with 



ui(A,A') = ln|A-A'| 
U2(A,A') = In 



( sinh"^ \/A ) - [smhr^ 



(17) 

(18) 
(19) 



We note that ui is translationally invariant in the variables A, A', but U2 is only invariant in the metallic regime 
(A <C 1). However, in the insulating regime (A » 1) U2 is negligible compared to ui. This fact suggests that we write 
U2 for the shifted variables 77 = A — A3, 77' = A' — A3 as 



Au2(r?,77'; A3) = In ■ 



where 



Then, the integral equation ( |16D is written as 
ds 



U2{v + + h) = u"{v,v') + Au2(77,7/; A3), 

(sinh-i VvTX^)^ - (sinh-i Vv^TX^)^ 



(sinh^^ ^)2 - (sinh"^ vV) 



7^2 



(20) 



(21) 



/ 

J — c 



ln| sinh(t - s)\ -|- In |i - s\ + -Au2{t,s) 



ospis) = 7i(A-|- A3,Ai) -)-7i(A + A3,A2) 
iT 



1 + A + A. 



+ Vi\ + \3), 



(22) 



where asp{s)d{s) = Psp{r] + X3)d{ri) with psp{X)d{X) = asp{x)dx, while sinh^t = 77 and sinh^ s = 77'; we have also 
extended the last integral to negative values since asp{s) = asp{—s). Finally, it turns out that Au2{t,s) can be 
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approximated by the product of the sum of two Lorentzians whose parameters are determined by the hmits s — > 0, 
s — » oo and s — X^. However, a posteriori numerical calculation showed that the contribution to the saddle point free 
energy Fsp{xi,X2, xy, <t{x)) coming from the Au2{t, s) term was negligible. So, in order to simphfy our calculations we 
neglect this term. Our last simplification is related to the interaction terms in asp{x): we consider only interactions 
between neighboring eigenvalues i.e., xi,X2 and X2,X3 and neglect interactions between xi and X3. As mentioned 
before, we will check our approximations by comparing < g > and var(f;) with those from the a model. 

Since the right hand side of Eq. ( p2| ) is linear in r and therefore asp{s), the saddle point free energy (Eq. (p^)) is 
quadratic in r, so Fgp can be written as 

Fsp =F° + {it)F' + ^-^F" (23) 
and the integral over r in Eq. ( p^ can be done exactly with the following result: 

nOO nOO pOG 

P{g) = / dxi dx2 / dx^e-^, (24) 

where 

^ F')Vi^°. (25) 



2F" 



Therefore once F^p is obtained from Ggpix), the calculation of the distribution Pijj) is reduced to a triple integration, 
Eq. (p^), which we compute numerically. 

Note that the above described method can in principle be implemented fully numerically, without making approxi- 
mations on AM2(i, s) or asp{x). The analytic approach allows us to identify the dominant terms and understand their 
change with disorder. The price we pay is that the expressions for the density and free energy are not valid for all 
possible values of Xi (0 < a;^ < cxi). In fact, our expressions for Gsp{x) and F^p are restricted to values of x-^ < 5. We 
will see, however, that it is enough to reach the region of main interest: the crossover regime. 



IV. RESULTS AND CONCLUSIONS 



We start comparing our calculations to the first two moments available for all degree of disorder from the a model 
[^3| . Figure |^ shows that the average and variance of g is now better in the crossover region, compared to the previous 
calculation of see Fig. |l|. On the other hand, the metallic regime is not described as good as in since we have 
neglected Au2 in Eq. ( |2C| ) and the interaction terms between Ai and A3. Actually, the approximations made in [l^ 
are quite good in this metallic regime, Fig. |l|, and we have included the results (squares) for the average and variance 
in Fig. H for completeness. Similarly, our restriction to values of X3 < 5 for asp{x) does not allow us to go throughout 
the insulating regime, since in this regime the conductance is dominated by the first two eigenvalues X2 ^ xi ^ 1. 
However, in our approach this region is very well described by two eigenvalues only as the density goes to zero. For 
completeness, we also include in Fig. || the results for such a calculation. Our main goal in this work is to develop 
approximations that are valid in the crossover region, and as shown in Fig. ^, our current approximations lead to 
very good results in the desired region. The complete Fig. ^ shows that the two limiting regimes can be described 
within the same formulation using different set of approximations. 

In figure ^ we show the evolution of P{g) as we change the disorder parameter F, evaluated within the current 
approximations valid in the crossover region. Suppose that we start decreasing the disorder in a sample i.e., we go 
from insulating to metallic behavior: at F = 0.5 we obtain a broad distribution with a bump at small values of g 
which is a reminiscence of the huge peak for << 1 in the insulating regime, see Fig. |^. Also, at 17 > 1, P{g) has 
an exponential decay as has been seen in numerical simulations [|[ 0, ^, |^ rather than the Gaussian cutoff at 5 > 1 
predicted in 10. As F is increased (F = 0.7), the bump at g << 1 disappears which makes < g > increase and at 
the same time, the decay at g > 1 becomes smoother . Finally, when we decrease the disorder even more (F = 1.0), 
we obtain a distribution which starts to look like a Gaussian distribution, as is expected in the metallic regime. 

In order to check our results in the crossover regime, we compare P{g) at F = 0.5 with the distribution obtained 
numerically by Plerou and Wang . Figure ^ shows a good agreement with the numerical simulation ^ , except 
at g=l where our P{g) has a peak. Others numerical simulations |0, ||, ^, |l^ do not show any cusp at g=l either. 
We believe that this behavior at g=l comes from our restriction X3 < 5 in the density asp{x), since as we go from 
the metallic to insulating regime large eigenvalues become more and more important. However, this peak eventually 
develops into the expected peak of the Gaussian in the metallic regime. On the other hand, an essential singularity 
of P{g) at g = 1 has been found in in the crossover regime, on the insulating side. So it is conceivable that the 
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peak present for F = 0.5 due to the existence of this singularity disappears, since the method developed here may not 
be valid in the presence of such singularities. Also, a non analytic behavior at g = 1 has been found numerically in 
quasi one and three dimensional systems in the crossover regime . 

In conclusion, we have developed a systematic method which allow us to calculate the distribution of the conductance 
in a quasi onc-dimensional system. In particular we were interested in the distribution of g across the crossover 
regime which had not been obtained before iQ . Separating out the two lowest eigenvalues and treating the rest as a 
continuum in the solution of the DMPK equation, we were able to obtain the evolution of P{g) as function the disorder 
parameter F. We have shown how P{g) develops from a broad flat distribution at the crossover regime to a Gaussian 
distribution in the metallic regime. Our results for the average and variance agree with the exact results from the 
non linear sigma model and the distribution P{g) in the crossover regime agrees well with the numerical calculation 
, including an exponential decay for g > 1 reported in numerical simulations. Together with the earlier results of 
how a log-normal distribution on the insulating side develops a cutoff at g = 1 and eventually becomes an asymmetric 
'one-sided' log-normal distribution near the crossover regime, this provides a relatively complete picture of how a 
log-normal distribution becomes a Gaussian via a highly asymmetric intermediate one-sided log-normal distributions 
when disorder is changed. The quantitative details at g = 1 at the crossover point remains an open question. 

Although our results are valid strictly only for quasi one-dimensional systems where the DMPK equation is valid 
(however, it has been shown that the restriction L ^ W can be relaxed (l8|), a qualitatively similar behavior of P{g) 
is found numerically in 2 and 3 dimensional systems in the crossover regime |l^. It is therefore possible that the 
method presented here could help understanding the generic behavior in higher dimensions at the crossover regime. 
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FIG. 1: (a) Mean and (b) variance of q calculated from the results in [g^ (squares) are compared with the correspondent 
results from the a model (solid line) [|l3|]. Horizontal arrows show schematically the different regimes of transport. Note that 
< g > and va.r{g) from |l4] are not calculated for the complete crossover region, since the deviation from the a model results 
becomes large. 




FIG. 2: Distribution P{g) in the metal (dashed line) and insulating (solid line) regime calculated by using the results in [ |l4| . 
As expected, a Gaussian distribution is obtained for the metallic case while, in a logarithmic scale, P{lng) follows a log-normal 
distribution for the insulating case. 




FIG. 3: Comparison of our results for the average and variance of the conductance, (a) and (b) respectively, with the non linear 
(T model (solid line) . Squares are the same as in Fig. 0. Dots are results obtained using the current approximations for the 
crossover region. Triangles are results obtained using two eigenvalues only. 
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FIG. 4: Evolution of the distribution of the conductance P{g) as the disorder parameter F is changed. Three cases are shown 
F = 0.5 (dots), F — 0.7 (squares) and F — 1.0 (diamonds). 



